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PREFACE 

This  report  describes  a method  for  estimating  inlet  velocities, 
discharge,  and  bay  levels  based  on  the  numerical  model  of  Seelig,  Harris, 
and  Herchenroder  (in  preparation,  1977).  This  method  for  predicting 
inlet  hydraulics  is  not  discussed  in  the  Shore  Protection  Manual  (SPM) 
(U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
1975).  The  work  was  carried  out  under  the  General  Investigation  of 
Tidal  Inlets  (GITI)  of  the  U.S.  Army  Coastal  Engineering  Research 
Center  (CERC) . 

The  report  was  prepared  by  William  N.  Seelig,  Research  Hydraulic 
Engineer,  under  the  general  supervision  of  Dr.  R.M.  Sorensen,  Chief, 
Coastal  Structures  Branch. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th 
Congress,  approved  31  July  1945,  as  supplemented  by  Public  Law  172, 
88th  Congress,  approved  7 November  1963. 


:V^L 


JOHN  H.  COUSINS 

Colonel,  Corps  of  Engineers 

Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


To  obtain 


To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C = (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K = (5/9)  (F  -32)  ♦ 273.15. 
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SYMBOLS  AND  DEFINITIONS 


A bay 
Ao 

Cl,  C2 

&bay 

4 max 
D 

g 

hfc 
hs 
Lb  ay 
Lin 

T > 

t 

0 


bay  surface  area  (square  feet) 

bay  surface  area  at  datum  (square  feet) 

coefficients  to  evaluate  Manning's  n (dimensionless) 

depth  of  bay  (feet) 

maximum  water  depth  in  inlet  (feet) 

Stillwater  depth  (feet) 

acceleration  of  gravity  (32.2  feet  per  second  squared) 

water  level  in  bay  (feet) 

water  level  in  sea  (feet) 

length  of  bay  (feet) 

length  of  inlet  (feet) 

forcing  wave  period  (seconds) 

time  step  used  in  model  (seconds) 

bay  surface  area  variation  parameter  (dimensionless) 
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A SIMPLE  COMPUTER  MODEL  FOR  EVALUATING  COASTAL  INLET  HYDRAULICS 

by 

William  N.  Seelig 

I.  INTRODUCTION 

This  report  describes  a method  for  estimating  coastal  inlet  veloc- 
ities, discharge,  and  bay  levels  using  the  simple  numerical  model  of 
Seelig,  Harris,  and  Herchenroder  (in  preparation,  19 77)  1 . The  model 
can  be  used  for  sea  level  fluctuations  caused  by  astronomical  tides, 
storm  surges,  seiches,  or  tsunamis.  A digital  computer  program  is  used 
because  of  the  large  number  of  computations.  A run  on  a CDC  6600  com- 
puter generally  costs  less  than  $5  for  a tidal  cycle. 

II.  PREDICTING  INLET  HYDRAULICS 
1 . Systems  Modeled  with  Computer  Program. 

An  inlet-bay  system  consists  of  a "sea"  (e.g.,  ocean  or  lake)  con- 
nected to  a "bay"  by  one  or  more  inlets  (Fig.  1).  The  computer  model 
will  predict  bay  levels,  inlet  velocities,  and  discharge  as  a function 
of  time  given  the  geometry  of  the  system  and  the  water  level  fluctua- 
tions in  the  sea.  It  is  assumed  ttiat  the  sea  is  much  larger  than  the 
inlet  and  bay  and  that  the  bay  is  large  compared  to  the  inlet. 

The  model  is  designed  for  systems  where  the  bay  water  level  rises 
and  falls  uniformly  throughout  the  bay.  This  occurs  when  the  wavelength 
in  the  bay  is  much  longer  than  the  longest  axis  of  the  bay: 


^^bay  ^ bay  > 

where 

Ip  = forcing  wave  period 

g = acceleration  of  gravity 

^bay  = depth  of  bay 
L bay  = length  of  bay 

2 . Procedures  for  Use  of  Computer  Program. 

Step  1.  Evaluate  the  inlet  geometry  by  using  maps,  charts,  hydro- 
graphic  surveys,  and  dredging  records  to  determine  the  depth  of  water 
throughout  the  inlet.  The  side  slope  of  the  inlet  at  mean  water  level 

^EELIG,  W.N.,  HARRIS,  D.L.,  and  HERCHENRODER,  B.E.,  "A  Spatially  , 
Integrated  Numerical  Model  of  Inlet  Hydraulics,"  GITI  Report  14,  U.S.  1 ' 
Army  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  Fort 
Belvoir,  Va. , and  U.S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Miss,  (in  preparation,  1977). 
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should  also  be  Measured.  Whenever  possible,  obtain  this  information  for 
the  time  of  interest  because  inlets  frequently  change  shape,  especially 
during  Major  storms. 

Step  2.  Construct  a flow  net  (series  of  cross  sections  and  channels) 
for  the  inlet  to  represent  the  Model  grid  (Fig.  1).  TTie  flow  net  and 
inlet  discharge  are  used  to  determine  bottom  friction  throughout  the  in- 
let. The  flow  net  is  drawn  by  approximating  the  average  path  (channel) 
that  water  follows  during  ebb  flow  and  floodflow.  Channel  boindaries 
are  drawn  along  these  paths  for  up  to  seven  channels.  A simile  inlet 
with  constant  depth  and  width  may  be  modeled  with  one  or  two  channels. 
Complex  inlets  require  approximately  three  to  seven  channels.  Channels 
should  have  the  smallest  spacing  in  deep  parts  of  the  inlet  where  flow 
will  be  highest.  Up  to  eight  cross  sections  should  then  be  drawn  per- 
pendicular to  the  channels.  The  first  cross  section  in  the  sea  and  the 
last  cross  section  in  the  bay  should  have  cross-sectional  areas  10  times 
larger  than  the  minimum  cross-sectional  area.  Cross  sections  should  be 
drawn  with  the  narrowest  spacing  near  the  minimum  cross-sectional  area 
section  where  friction  in  the  inlet  will  be  high. 

Step  3.  Measure  the  surface  area  of  the  bay  at  the  mean  water  level, 
Ao,  from  charts  or  aerial  photos.  For  most  bays  the  surface  area  changes 
as  the  bay  water  level  rises  and  falls  because  sections  are  flooded  at 
high  water  levels.  If  the  bay  area  change  is  significant,  a bay  area 
variation  parameter,  3,  is  used  to  account  for  area  of  the  bay,  A^^, 
at  any  water  level  in  the  bay,  h£,  using  the  relation: 

N)ay  = Ao^  + » (2) 

where  Aq  is  the  bay  surface  area  at  datum,  usually  mean  low  water  (MLW) , 
mean  sea  level  (M5L)  , or  mean  water  level  (ML) . 


Step  4.  Specify  the  seawater  level  fluctuation  as  a function  of  time 
for  the  period  of  interest.  Tide  tables  will  give  an  estimate  of  the 
astronomical  tide.  Water  levels  can  also  be  measured  by  a tide  gage  and 
stilling  well  (Seelig,  1977)  2.  Corps  of  Engineers  and  National  Oceanic 
and  Atmospheric  Administration  (NOAA)  gages  located  at  numerous  points 
along  the  coast  may  also  provide  the  desired  water  level  information. 

In  this  co^mter  program  either  the  tide  may  be  expressed  as  a sinusoidal 
wave  with  a period  and  a^ilitude  or  the  levels  may  be  described  by  instan- 
taneous sea  level  measurements  at  a constant  sailing  rate. 

Step  5.  Determine  the  time  step  of  input  to  the  model  for  use  in 
computations.  As  a lower  limit,  the  time  step.  At,  should  be: 


At 


(3) 


SEELIG,  W.N. , "Stilling  Well  Design  for  Accurate  Water  Level  Measure- 
ment," TP  77-2,  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  Fort  Belvoir,  Va. , Jan.  1977. 


where  L^n  is  the  length  of  the  inlet  and  ^rriax  is  the  maximum  water 
depth  in  the  inlet.  A longer  time  step  can  be  used  for  most  tidal  in- 
lets, and  as  an  upper  limit,  the  time  step  should  be  one-hundredth  of 
the  forcing  wave  period. 

Step  6.  Document  all  input  data  using  the  computer  format  shown  in 
the  appendix.  As  a first  estimate,  set  the  flood  and  ebb  entrance  and 
exit  loss  coefficients  to  equal  one  (CDF  = 1.0  and  CDE  = 1.0).  As  a 
first  approximation.  Manning's  n can  be  evaluated  by  the  relation: 

n = Cl  - C2  D , (4) 

where  D is  the  local  inlet  Stillwater  depth.  For  depths  greater  than 
4 feet  and  less  than  30  feet.  Cl  = 0.03777  and  C2  = 0.000667;  for  depths 
less  than  4 feet.  Cl  - 0.0550  and  C2  = 0.005.  The  n for  each  grid  may 
be  different  if  C2  i 0. 

Step  7.  For  use  with  periodic  forcing,  run  the  program  for  several 
sinusoidal  cycles  having  the  period  and  amplitude  of  the  long  wave  of 
interest  to  approximate  the  hydraulic  characteristics  of  the  inlet-bay 
system.  A sinusoidal  tide  is  specified  in  the  model  by  giving  the  forc- 
ing period,  T,  in  hours  and  the  wave  amplitude,  \Q,  in  feet,  on  card 
type  3 and  by  setting  NPTS  = 0 on  card  type  8 of  input  to  the  program 
INLET.  Set  ITABLE  = 1 to  obtain  tables  of  instantaneous  hydraulics  at 
points  throughout  the  water  level  cycle  and  set  IPLOT  = 1 to  obtain  a 
plot  of  predicted  inlet  velocities  and  discharge  at  sequential  bay  levels 
These  outputs  will  indicate  the  importance  of  the  terms  in  the  equation 
of  motion  describing  water  motion  in  the  inlet.  If  temporal  acceleration 
is  small  during  most  of  the  water  level  cycle,  then  startup  transients 
will  be  small  and  the  first  or  second  cycle  will  contain  little  transient 
effect  (NCYCLES  = 1 or  2 in  input  data) . However,  if  temporal  accelera- 
tion is  significant  during  more  than  25  percent  of  the  cycle,  approxi- 
mately four  cycles  of  model  operation  are  required  to  eliminate  startup 
transient  effects  (NCYCLES  = 4).  For  aperiodic  use  such  as  with  storm 
surges  or  rapidly  varying  wave  size  (e.g.,  tsunamis),  run  the  model  for 
the  water  level  for  approximately  10  hours  before  the  time  of  interest 
to  build  up  initial  conditions  in  the  model  similar  to  the  prototype. 

Step  8.  Calibrate  the  computer  model  by  varying  Manning's  n or 
flood-  and  ebb- loss  coefficients.  The  seawater  level  fluctuation  can  be 
specified  as  a sinusoidal  wave  or  in  terms  of  an  equal  time  series.  For 
an  equal  time  series,  start  and  stop  the  series  when  the  seawater  level 
is  at  zero  so  that  one  or  more  complete  cycles  are  described.  Use  at 
least  20  points  to  describe  each  cycle.  The  sampling  interval  in  minutes 
TDEL,  and  the  number  of  points,  NPTS,  must  be  specified  on  card  type  8 
and  the  water  level  data  on  card  type  9. 

The  model  is  calibrated  using  short  periods  of  field  observations  by 
first  comparing  observed  and  predicted  mean  water  velocities,  if  avail- 
able, at  the  minimum  cross-sectional  area  region  of  the  inlet.  If  the 
predicted  velocities  are  higher  or  lower  than  observed,  then  the  value 
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of  n can  be  increased  or  decreased  accordingly.  When  the  computer 
model  has  been  satisfactorily  calibrated  to  predict  inlet  velocities, 
predicted  bay  water  levels  should  be  checked  against  measurements  to 
assure  that  levels  are  being  modeled  correctly.  If  inlet  velocities 
are  not  available,  bay  levels  can  be  used  to  calibrate  the  model. 

Step  9.  If  additional  prototype  data  are  available,  these  data 
should  be  used  to  verify  that  the  model  adequately  predicts  inlet  and 
bay  hydraulics. 

Step  10.  At  this  point  the  computer  program  is  ready  to  use  for 
prediction.  Examples  of  the  use  of  the  computer  program  are  presented 
in  the  following  section.  Input  and  output  data,  and  computations  are 
in  U.S.  Customary  units. 

III.  EXAMPLES  OF  COMPUTER  PROGRAM  PREDICTION 
1 . Cabin  Point  Creek,  Virginia. 

Cabin  Point  Creek  is  a shallow  natural  tidal  inlet  that  connects  a 
bay  to  the  lower  Potomac  River  (Fig.  2)  where  the  mean  tidal  range  is 
approximately  1.5  feet. 

In  this  example,  the  model  was  calibrated  with  prototype  river  and 
bay  levels  and  the  calibrated  model  was  then  used  to  predict  inlet  veloc- 
ities, discharge,  and  bay  level  for  a second  inlet  added  to  the  system. 
The  procedures  for  using  the  model  are: 

(a)  The  inlet  cross  section  was  measured  (Fig.  3)  on  24 
May  1976,  and  is  assumed  to  be  representative  of  the  1,900- 
foot-long  inlet. 

(b)  The  inlet  is  modeled  using  a grid  system  of  three 
channels  and  two  identical  cross  sections  (Fig.  3)  at  either 
end  of  the  inlet. 

(c)  The  bay  area,  Ac,  measured  from  a T^-minute  U.S. 

Geological  Survey  (USGS)  topographic  map,  was  3.5  * 106 
square  feet.  For  an  increase  in  bay  water  elevation  of 
0.25  foot,  the  bay  surface  area  increases  approximately  5 
percent  because  of  marsh  flooding.  The  bay  area  variation 
parameter,  &,  can  be  determined  from  this  information 
using  equation  (2),  rearranged  as: 
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or,  in  this  case, 

«-51s'us  - » ■ »•’ 


Oepth  (ft) 


Bay  Wafer  Level 
Measurements 


Cobin 

X^Poinf 

^ VIRGINIA 

I 

Design  Inlet 
(inlet  2) 


POTOMAC 

RIVER 


Naturol  Inlet 
( inlet  I ) 


River  Level 
Measurements 


Approx.  Scale  (m) 

I — i l — l =3 

SOO  1.000 

Approx.  Scole  (ft) 


(d)  River  water  levels  were  measured  at  30-minute  intervals 
using  a stilling  well  located  near  the  inlet  mouth  (Fig.  2). 


(e)  The  time  step  was  estimated  as: 

*«-  1900  . 

✓32.2 x 2 

(f)  Loss  coefficients  were  specified  as  CDF  = CDE  = 1.0,  and 
Manning's  n was  estimated  as  n = 0.055  - 0.005  D (recommended  for 
depths  less  than  4 feet) . 

(g)  A preliminary  computer  run  using  a sinusoidal  river  tide 
showed  that  the  inlet  is  controlled  by  friction  effects  and  that 
tenoral  acceleration  is  not  iuportant. 

(h)  The  model  was  then  run  using  the  measured  river  water 
levels  to  force  the  model  (Fig.  4).  It  was  determined  that  the 
model  adequately  predicted  bay  levels. 

(i)  No  additional  prototype  data  are  available  for  verifi- 
cation of  the  model. 

(j)  The  model  is  now  available  to  use  for  predictions  of  in- 
let hydraulics.  In  this  example,  a second  inlet  (inlet  2)  , is 
being  considered  for  this  site,  so  the  model  is  used  to  predict 
hydraulics  for  the  system  with  two  inlets  (Fig.  2).  Procedures 
(a)  and  (b)  are  repeated  for  the  second  inlet.  In  this  case,  the 
second  inlet  is  modeled  by  one  channel  and  two  cross  sections  so 
that  the  inlet  has  a length  of  300  feet,  a width  of  50  feet,  and 
a depth  of  4 feet.  These  inlet  data  are  put  into  the  conputer 
format,  added  to  the  program  deck  for  the  natural  inlet,  and  re- 
run to  predict  conditions  for  the  proposed  two-inlet  system.  The 
numerical  model  predicts  that  addition  of  the  second  inlet  would 
increase  the  tidal  range  and  the  tidal  prism  in  the  bay  and  would 
cause  water  velocities  in  inlet  1 to  decrease  (see  Table). 


Table.  Predicted  Cabin  Point  Creek  hydraulics. 


Tide 

24  and  25 
May  1976 

Model  prediction 
for  second  inlet 

Inlet  1 

Inlet  1 

Inlet  21 

Bay 

(range  in  ft) 

0.36 

1.49 

Hi 

Ebb 

(maximum  velocity  in  ft/s) 

-0.6 

-0.3 

Flood 

(maximum  velocity  in  ft/s) 

0.9 

0.3 

1L  = 300  feet,  B = 50  feet,  D = 4 feet. 


NOTE:  Tidal  range  in  the  sea  is  1.49  feet. 
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Level  (ft) 
arbitrory  datum 


24  May  Time(hr)  25  May 

Figure  4.  Cabin  Point  Creek  sea  and  bay  levels. 
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2.  Pentwater  Inlet,  Michigan. 

Pentwater  Inlet  is  an  example  of  a Great  Lakes  inlet  controlled  by 
vertical-walled  jetties  along  the  entire  2,000-foot  channel  (Fig.  5). 
Meteorologically  generated  seiches  of  Lake  Michigan  are  the  primary  water 
level  fluctuations  causing  reversing  currents  in  the  inlet.  A model  of 
Pentwater  will  be  calibrated  and  used  to  estimate  hydraulic  response  of 
the  inlet  to  simultaneous  lake  seiching  and  river  inflow.  The  procedures 
used  in  this  modeling  are: 

(a)  A hydrographic  survey  of  the  inlet  is  used  to  describe 
the  inlet  geometry. 

(b)  The  inlet  is  modeled  using  one  channel  and  six  cross 
sections . 

(c)  The  bay  surface  area,  measured  from  a hydrographic 
chart,  is  1.81  x 107  square  feet.  The  bay  area  does  not 
change  with  bay  water  level  because  the  bay  has  steep-sided 
slopes,  so  8 = 0. 

(d)  Lake  Michigan  water  level  measurements  used  to  force 
the  model  were  taken  at  5-minute  intervals  on  a tower  located 
adjacent  to  Pentwater  Inlet. 

(e)  The  model  time  step  used  is: 

At  = 2000==^  _ go  seconds 

✓32.2  x 15 

(f)  Loss  coefficients  were  specified  as  CDE  = CDF  = 1.0, 
and  Manning's  n was  estimated  by  n = 0.03777  - 0.000667  D 
(recommended  for  depths  greater  than  4 feet  and  less  than  30 
feet) . 

(g)  A preliminary  run  showed  that  temporal  acceleration  is 
an  important  term  in  the  inlet  equation  of  motion  for  Pentwater 
Inlet  (Fig.  6).  Therefore,  several  forcing  cycles  of  model 
operation  before  the  time  of  interest  are  necessary  to  eliminate 
transient  terms  due  to  startup  conditions. 

(h)  The  model  is  calibrated  by  using  Lake  Michigan  levels 
to  force  the  model.  An  initial  run  showed  that  predicted  bay 
level  fluctuations  adequately  modeled  observed  levels  (Fig.  7). 

(i)  The  model  was  not  verified. 

(j)  The  model  was  used  to  predict  inlet  velocities,  dis- 
charge, and  bay  levels  for  a 2-hour  forcing  wave  with  an 
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Importance  of  Terms  in  the  Eqoation  of  Motion 
( normolized  by  the  lorgest  term ) 


Bolton  Stress 

Pressure  Head 
Convective  Acceleration 
Tenporal  Acceleration 


Tronsient  Effects 


Monochronatic 


Figure  0.  Pentwater  Inlet  model  prediction  of  monochromatic  forcing 
(for  a 2-hour  have  with  a 9.1-foot  amplitude). 
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Water  Level  (ft) 


Measured  Offshore  Levels 


Time  (hr) 

Figure  7.  Pentwater  Inlet  model  calibration. 


amplitude  of  0.10  foot  and  a discharge  into  Pentwater  Lake  of 
2,800  cubic  feet  per  second  from  the  Pentwater  River.  The  model 
predicted  an  average  bay  water  surface  elevation  of  0.13  foot 
higher  than  the  mean  lake  level,  a bay  water  level  fluctuation 
range  of  0.25  foot,  and  a prism  of  water  of  4.6  * 106  cubic  feet 
caused  by  the  seiche  (Fig.  8).  The  inlet  would  always  be  in  ebb 
flow  due  to  river  influence  with  a maximum  velocity  of  -2.7  feet 
per  second  and  a minimum  velocity  of  -0.1  foot  per  second.  Head, 
friction,  and  temporal  and  convective  acceleration  are  important 
in  the  inlet  equation  of  motion. 

IV.  SUMMARY 

A computer  program  (INLET)  based  on  a numerical  model  (Seelig,  Harris, 
and  Herchenroder,  in  preparation,  1977) 1 is  presented  for  prediction  of 
hydraulics  where  one  or  more  inlets  connect  a bay  to  a sea.  Two  examples 
are  given:  (a)  A tidal  inlet  forced  by  an  astronomical  tide  where  inlet 
channel  friction  is  the  dominant  term  in  the  equation  of  motion;  and  (b) 
a Great  Lakes  inlet  with  river  inflow  forced  by  lake  seiching  where  head, 
friction,  and  temporal  and  convective  accelerations  are  important  at 
different  points  in  the  water  level  fluctuation  cycle.  The  model  can 
also  be  used  for  forcing  other  water  level  fluctuations,  such  as  from 
storm  surges  or  tsunamis. 

Another  computer  program  (INLET2)  is  available  for  more  complex 
systems  of  interconnected  inlets,  bays,  and  seas.  INLET2  is  an  expanded 
version  of  INLET.  Documentation  and  conputer  card  decks  for  INLET2  are 
available  from  the  Automatic  Data  Processing  Division  (CERDP) , Coastal 
Engineering  Research  Center  (CERC) . 

Details  on  model  development  and  application,  including  additional 
examples,  are  reported  by  Seelig,  Harris,  and  Herchenroder  (in  prepara- 
tion, 1977) 1 . 


Seelig,  Harris,  and  herchenroder,  op.  cit.,  p.  7. 
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Pressure  Head 

Convective  Acceleration 

Temporal  Acceleration 


Time  (hr) 


Figure  8.  Predicted  Pentwater  Inlet  velocities,  discharge,  bay 

levels,  and  relative  magnitude  of  terms  in  the  equation 
of  notion. 


APPENDIX 


COMPUTER  PROGRAM  DOCUMENT  AT I OX  (INLET) 


1.  Program  Description. 


The  numerical  model  to  predict  inlet  hydraulics  is  programed  in 
FORTRAN  for  a CDC  6600  computer.  The  simultaneous  differential  equations 
are  solved  by  a variable  time  step  Runge-Kutta-Gill  marching  procedure. 
The  organization  of  the  computer  program  is  shown  in  Figure  A-l.  A brief 
description  of  each  routine  follows: 


INLET  is  the  main  routine  which  controls  input-output  and  calls 
subroutines  to  execute  a specific  task.  Figure  A-l  summarizes  control 
throughout  the  program.  The  program  is  organized  to  accept  to  three 
inlets  connecting  the  bay  to  the  sea,  up  to  seven  channels  for  each  in- 
let, and  up  to  ei  ght  cross  sections  (seven  grids  long) . 


Subroutine  HELM  uses  an  iterative  method  of  estimating  the  natural 

for  the  inlet-bay  system  by 


punping  period  or  Helnholtz  period,  T 


a 


neglecting  friction  in  the  inlet  to  give: 


(L in  + L')  A 


bay 


&V 


where  L*  is  added  inlet  length  due  to  radiation,  and  where  L’  is 
given  by: 


Subroutine  RKGS  is  a routine  to  solve  simultaneous  differential 
equations.  This  sii»routine  was  adapted  from  the  scientific  subroutine 
package . 

Subroutine  SETEQ  evaluates  the  right-hand  side  of  the  equation  of 
motion,  one  for  each  inlet,  and  the  continuity  equation  between  the  inlet 
and  bay  for  each  step.  This  routine  also  evaluates  the  relative  rank  of 
the  four  terms  in  the  equation  of  motion  for  flow  in  each  inlet. 

Sii>routine  LEVEL  determines  the  water  level  in  the  grids  at  each 
time  step.  The  routine  interpolates  the  level  between  the  sea  and  bay 
based  on  the  relative  amount  of  friction  in  each  grid  cell. 

Sii) routine  TPWRTE  writes  hydraulic  results  from  each  time  step  on  a 
tape  or  disc,  so  that  this  information  can  be  used  later  by  the  output 
routines. 

Subroutine  TABLE  outputs  a table  of  instantaneous  hydraulics  each 
time  the  routine  is  called. 
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Subroutine  SF.A  determines  the  water  level  in  the  sea  as  a function 
of  time  either  for  a given  sine  wave  or  by  interpolating  equal  time- 
series  data. 

Subroutine  WT1  determines  the  grid-weighting  function  by  assuming 
that  the  flow  is  distributed  across  each  section  so  that  friction  is 
minimized.  This  routine  allows  flow  to  cross  channel  boundaries,  but 
assumes  that  this  flow  will  be  small,  so  the  flow  is  neglected  in  the 
equation  of  motion.  This  weighting  function  is  recommended  for  general 
use. 


Subroutine  WT2  is  similar  to  WT1,  except  that  flow  is  not  allowed  to 
cross  channel  boundaries  and  that  flow  is  distributed  in  each  channel  so 
that  friction  is  minimized. 

Subroutine  1VT3  determines  the  weighting  function  so  that  flow  is 
distributed  equally  in  all  grids.  This  is  generally  unrealistic,  since 
it  will  be  difficult  to  visually  draw  this  grid  system.  However,  this 
routine  is  useful  since  it  provides  an  upper  limit  on  frictional  effects 
and  therefore  gives  a lower  limit  of  bay  levels  and  inlet  velocities. 

This  weighting  can  be  used  to  model  simple  geometry  inlets  where  only 
one  channel  is  used  to  represent  the  inlet. 

Subroutine  CRIT  prints  a table  of  critical  instantaneous  hydraulics 
(i.e.,  at  high  water,  low  water,  maximum  velocity,  and  maximum  discharge). 
This  table  is  determined  by  storing  a summary  of  conditions  for  each  time 
step,  then  scanning  this  list  for  critical  values. 

Subroutine  GRPHC  plots  mean  inlet  hydraulics  by  scaling  hydraulics 
in  storage  and  plotting  the  time  interval  requested  on  a digital  x-y  pen 
plotter. 

Subroutine  READ IN  is  used  by  GRPHC  to  read  data  in  storage  and  scale 
values  for  plotting. 

2.  Program  Input. 

The  computer  program  (INLET)  requires  the  following  input  of  one  deck 


for  each  inlet-bay  system: 

Card 

type 

Variables 

Format 

Description 

1 

ALABL1 

4A10 

first 

line  of  title 

ALABL2 

4A10 

second 

line  of  title 

2 

N INLET 

5110,  2F10.5,  110 

number 

of  inlets 

NCYCLES 

I PLOT 

number 

I PLOT  = 

of  cycles 

= 1 for  plot  of  results 
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Format 


Description 


Card  Variables 
type 

HIT 


I tabij: 
Cl,  C2 


I COW 

5 3F10.5,  E10.4, 

3F10.5,  2F5.1 

T 

DELT 

AO 

AB 

BETA 

IETA 

OINFLO 

CDF 

CDE 

4 2110,  F10.0 
IC 

IS 

QIJiT 

5 (one  card  10X,  7F10.5 

per 

section) 

A' 


6 (one  card  10X,  7F10.5 
per 

section) 

B' 


weighting  type 

IKT  = 1 flow  distributed  to  mini- 
mize (1  in  card  col.  40) 

ITABLE  = 1 for  tables  of  instanta- 
neous hydraulics 

Manning's  n evaluated  by: 
n = Cl  - C2  * D;  where  D is  still- 
water  depth.  If  blank  default 
values  of  Cl  = 0.03777  and 
C2  = 0.000667  are  assumed. 

I CONY  = 1 (1  in  card  col.  80) 


forcing  period  (hours) 
approximate  time  increment 
forcing  wave  amplitude  (feet) 
bay  area  at  datum  (square  feet) 
bay  area  variation  parameter 
inlet  side  slope  D(z)/D(y) 

bay  inflow  from  sources  other  than 
the  inlet  (cid>ic  feet  per  second) 

an  e^>irical  flood-loss  coefficient 
an  eiqiirical  ebb-loss  coefficient 

nud>cr  of  channels 

niaber  of  cross  sections 

estimated  inlet  discharge  at  the 
time  the  model  starts 


cell  cross-sectional  areas  at  the 

ends  of  each  cell  at  datum  (square 

feet)  (see  Fig.  A-2)  1 


grid  cell  widths  for  the  end  of 
each  cell  (feet)  (see  Fig.  A-2) 
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input  data. 


Format 


Description 


Card  Variables 

type  

7  (one  less  10X,  7F10.5 
card  than 
sections) 

L'  lengths  of  the  sides  of  cells  (see 

Fig.  A-2)  (one  less  card  than  num- 
ber of  sections;  one  more  value  per 
card  than  the  number  of  channels) 


For  card  types  5 to  7,  there  will  be  one  card  for  each  cross  section 
of  the  inlet.  The  first  card  will  be  for  the  first  cross  section,  i.e., 
the  section  closest  to  the  sea,  and  the  last  section  is  adjacent  to  the 
bay.  The  first  value  on  each  card  will  correspond  to  the  first  channel 
adjacent  to  land;  the  last  value  on  each  card  will  correspond  to  the 
last  channel  also  adjacent  to  land  (Fig.  A-2). 

For  more  than  one  inlet  connecting  the  bay  to  the  sea,  repeat  card 
types  3 to  7 for  each  additional  inlet. 

Card  Variables  Format  Description 

type  

8  TDFL  34X,  F6.2  water  level  sampling  interval 

(minute) 

NPTS  6X,  13  number  of  sample  points  = 0 for  no 

data 


9  (optional  -no  cards  if  NPTS  = 0 from  card  type  8) 

Y eight  water  level  values  per  card, 

as  many  cards  to  include  NPTS  points; 
start  the  model  at  a time  when  the 
sea  level  is  zero.  Use  25  or  more 
points  per  forcing  cycle  for  best 
results;  i.e.,  levels  at  30-  or  15- 
minute  intervals  for  a 12-hour  tide. 
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(optional--two  plot  cards,  first  card  used  only  if 
I PLOT  = 1 on  card  type  1) 


XO 

XF 

SCALX 

YLO 

YL 

YLSCAL 


8F10 . 5 ,/ , 3F10 . 5 , 110 

starting  time  of  plot  (hours) 
ending  time  of  plot  (hours) 
time  scale  (hours  per  inch) 
minimum  value  of  water  levels  (feet) 
overall  height  of  plot  (inches) 

scale  of  water  level  height  (feet 
per  inch) 
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Format 


Description 


A 


Card  Variables 

type 

YRO  minimum  flows  (thousand  cubic  feet 

per  second) 

YRSCAL  scale  of  flows  (thousand  cubic  feet 

per  second  per  inch) 

Second  card 

YVO  minimum  velocity  (feet  per  second) 

YVSCAL  scale  of  velocities  (feet  per 

second  per  inch) 

SCALE  scale  factor  for  total  plot  size 

IQ  IQ  = 0 for  no  plot  of  inlet 

discharge 

11  If  a plot  is  requested,  repeat  card  types  8 and  9 for  observed  bay 
levels  to  compare  with  predictions  (card  type  8 required;  use 
NPTS  = 0 for  no  observed  bay  levels) . Only  one  set  of  card  types 
10  and  11  will  be  required  for  plotting  even  though  the  system 
modeled  may  have  more  than  one  inlet. 

12  End  of  file  card. 


The  inlet  data  for  a computer  run  of  Masonboro  Inlet,  North  Carolina, 
are  shown  in  Figure  A- 3. 


Program  Output. 


The  types  of  output  include:  (a)  A summary  table  of  grid  dimensions, 
input  parameters,  and  the  Helmholtz  period  of  the  system  estimated  assum- 
ing there  is  no  friction  in  the  inlet;  (b)  (optional)  summary  tables  of 
instantaneous  inlet  hydraulics;  (c)  (optional)  a pen  plot  of  inlet  hydrau- 
lics; and  (d)  a table  summarizing  critical  points  throughout  model  opera- 
tion, such  as  high  water,  low  water,  point  of  maximum  discharge,  and 
maximum  velocity.  Samples  of  input  and  output  for  the  Masonboro  Inlet 
run  are  given  in  Figures  A-4,  A-5,  and  A-6. 


suter  Program. 


A listing  of  the  computer  program  (INLET)  follows  the  sample  output. 

The  program  was  written  in  FORTRAN  IV  for  a CDC  6600  computer  with  plotter. 
Control  cards,  plotting  instructions,  and  file  controls  may  have  to  be 
changed  for  other  computers.  If  no  plotter  is  available,  the  subroutine 
GRPHC  and  the  call  to  the  subroutine  in  the  main  program  may  be  removed. 
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Figure  A-3.  Sample  of  input  data  for  a computer  run  of 
Masonboro  Inlet,  Nort'i  Carolina. 
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-.«* 

-.*9 
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.1* 

.•2 

•.2* 

1.T0 

2.A0 

2.33 

2.«0 

2.50 

2.41 
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.5" 

ft. .'ft 

-.5f 
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• 1 .*2 

• 1.00 

-l.A* 

•••03 

•.ft* 
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• 3* 

.•3 
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• .* » 

2.'- 

7.  3t 

>.JN 

>.4* 

«.*7 

1.5* 

1.10 

.** 

• 10 

-.•ft 

-.•»  -1.3" 

Figure  A-4.  Sample  output  from  INLET  (summary  table  for 
Masonboro  Inlet  input  data). 
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1177?. 

1 1 77?. 
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1 177?. 
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.2| 

• ? 1 
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.21 

.21 
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EPIC 

• no 
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."2 

.10 
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.08 

if** 

»CC« 

.6  C0*«v  *cc * 

12. « 

HE»"«  -106 

.0  E a I C » 

67.0 

“E  »"•  VE|  OCITV 

AT 

Up  i»f» 

SEC’IO1  * 

2.97  77/StC  »“7N» 

18  2*. 71  E 72 

Figure  A-5. 


Sample  output  from  INLET  (summary  table  of  instantaneous 
hydraulics  for  Masonboro  after  6 hours  of  model  time). 
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11.13# 

-.279 

o#o*o 

,5«2 

-3.429* 

.59  . 1 it) 
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. a*5 

-3.433* 
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1 1.566 
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.52.606 
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.100 
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1 1 ."»« 
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.0«7 

-3.420* 
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-.°33 

0.0*0 

-.*)a 
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•1  , a95 
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• • 0 7 i 

-.923 
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-.•1? 

0.0*0 

•1.295 

1 .00** 

55.999 
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1.153 

0.000 

.105 

2.994* 

50.979 

17.160 

1 .25* 

0.0*0 

. 2H3 

3.020* 

52.006 

1 7, SO* 

1 .354 

0.000 

.30? 

3 • 0 3b* 

52.885 

17,687 
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17.77* 
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53.685 
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1 .*95 
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.*7  2 

3.0*4 
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iT.UHO 

1 .*.40 
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.719 
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53.719 

18,056 
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0.0*0 
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2. 99a 
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18.111 
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55.968 
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,99a 
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19.13# 
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1.003 
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21.779 
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.52,628* 

21.990 
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-2.921* 
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22.000 
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.?*« 

0.0*0 

1 .0*4 

-3.415* 

.56. #76 

21.000 

.15* 

0.000 

.970 

-3.4?9* 

.56.18# 

21.111 

. naa 

o.ooo 

,H7b 

-3.440* 

.55.836 

21.221 

-.0*7  • 

0.000 

.7*0 

-3.449* 

-5S.U60 

21.11# 

-.17* 

0.0*0 

.004 

-3.45** 

-55.0## 

21.«#5 

-.2*9 

0.000 

.507 

•3.469* 

-5#. 588 

21.956 

-.400 

0*000 

.409 

-3.401 • 

•5U.092 

21.667 

-.513 

o.ooo 

.390 

•3.401* 

-55.57# 

21.779 

- • *2* 

o.ooo 

.2q0 

-3.463* 

.55.063 

21.990 

-.7«1 

0.000 

.109 

-3.462* 

-52.516 

2#. 000 

-,•449 

0.000 

.007 

•3.45a* 

.51.870 

2#. Ill 

-.951 

0.000 

-.015 

-3.435* 

•51 .063 

2# . 221 

•1.052 

0.000 

-.117 

•3.409* 

-50.167 

25.000 

• 1.390* 

0.0*0 

-.055 

-2.599 

•A5.9#« 

• CRITICAL  POINT  VALUt 

Figure  A-6.  Sample  output  from  INLOT  (table  of  critical  points  for 
the  model  time:  high  water,  low  water,  etc.,  for 
Masonboro  Inlet). 
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1 97 1 

FOH«AT<S».  (niP7»(F7)  l,l».ll)F10.*l 

Inlet 

1 34 

1979 

F0M«»T(4».(LfN(FT)  (.2l.IOP16.il 

INLtT 

139 

1975 

Fl'P-AT  (5*  , (6  C.  IOF.  IflF  10.61 

INLtT 

US 

10 

CONTINUE 

INLtT 

136 

C eind  **m  *"0  -IoTh  *t  t*f  «ini><u«  section 

INLtT 

137 

A*IwJ(NI )*99,e«<2 

INLtT 

136 

DO  109  !■ 1 » 1 6 

inlet 

139 

AA«n, 

INLtT 

190 

B«*0# 

INLtT 

191 

DO  10*  .1*1.  IC 

INLtT 

192 

AA*AA*M I. J) 

INLtT 

193 

10* 

HH*RH  + H f I • J ) 

INLtT 

194 

IMU.Ct,»hM(n1)j  60  TO  109 

INLET 

1 9 S 

( '*  I ) * A A 

INLET 

196 

BMI*I  (*I  )*H6 

inlet 

197 

loo 

COM  J*UF 

INLtT 

196 

1110 

COM  J*uF 

INLET 

199 

C tST A Tf  T Hf  IHtT-M*  Pt»10D 

INLET 

ISO 

C*L|.  Mtl«fTHtLM.»H,CU»L) 

INLtT 

151 

7»7F»T/THt|_n 

INLtT 

1 S2 

»4lTt(6.*0l)  T,TnFl“.TF-TF 

INLtT 

1 S 3 

*01 

FO»x*T(t«,  (FciUCInS  pERI()04  l,FT,*<  I NUURSl. 

INLtT 

154 

!/•!*«  llHtt.X(0*PI)OI)B  |»f  *,*•  t HOURS  1«  / 

INLtT 

155 

1 1>.  PF/T"!  I.IOK.P*,*) 

INLtT 

156 

► 417  f,  (6.1557)  < ( J.LFNGTM(J) .LORLl J)) • J»1 .NINLET) 

INLtT 

157 

1 337 

F06"*T(  .|*.  (INLtT  LENGTH  FOOtO  LtNGTNl.  (/,«*.  I *.  1 « . 

INLtT 

156 

1 76. 1 1 )) 

INLtT 

159 

74745606, 

INLtT 

160 

CFLL  “AGS(t  N0.0lLT,NlNLET.(itNFL0.1TF8LE.T) 

Inlet 

161 

OELT4tND/FLOFT(NT> 

Inlet 

162 

DO  fib*  Nj4l.Nli.Lt7 

INLtT 

163 

NM4MS 

I Nl.  £ T 

164 

«»lTt(6.**6*)  Nt 

INLET 

165 

*206 

F OH"fT{//,  1 OX.  ISU"«I»V  TFOLF  OF  HVD4FULICS  INLE71.H1 

INLtT 

166 

IUN77»N1«6 

INLET 

167 

C»ll  C4|7(N7,DtL7.  !')NT7.7.NC7CLIS) 

Inlet 

166 

IFtrpLOT.tO.l.FNO.Nt.FO.l)  C*LL  PL07SC4UF.iOOO.51 

inlet 

169 

1F1 JPL07.E0.1)  CALL  G4PHC(*L*HL1.»LF«L*.0EL7.!UNJ7,NJ) 

INLt  T 

1 7C 

lE(TpLO7,tO.l,»N0.Nl,fO.NINLE7)  CALL  PL07 ( 0 . . 0 . . 444) 

INLtT 

171 

**64 

CON7  JNUF 

INLtT 

172 

5 704 

INLtT 

174 

END 

INLtT 

179 

34 


BEST  AVAILABLE  COPY 


SUbRUUTTNF  »KGS(END,DFLT,NlNLFT,UINFLO.ITAaLE,T) 

C ROUTINE  TO  SOL  VF  A SET  OF  5I"OTanEOUS  DIFFERENTIAL  FOUATIONS 
C AOARTFD  FROM  ScTFNTTFIC  SUBROUTINE  PACKAGE.  IBM,  1R70 
CrjMMUs/NUM|/Y(5),ntpv(s),x,MT,IwT,ZtT*,HS 
COMhON/RIIMU/RMk(  J,<i) 

01  "FUSION  AIIX(e,S),»(S),B(8},Cta),PRMTt5),AMlMI(l) 
NOI-taNlNLF  T»1 
RRMT(l)  • 1. 

RRMT(2)«tMO 
PR'*T  ( UbOELT 
PRMT(R)  » ,1 

I'lT.C.T. 56000.)  0ELTH356U0. 

IFtT.LF.jMOOO.)  0FlTBpT/». 
no  tl?<!  JNal.NlNLET 
vf jm)po.oi 
1112  OF  Rv ( Jn) bo ,no 1 

PFm»(NUIm).i ,o. FLOAT (MINUET )»0, 001 
ro  t i»i,noim 

1 All*  ft*.  I >30,06666666T30ERY  ( I) 

XbPrmT ( ) ) 

xtNo«PwMT (2) 

HBPR“T ( 3) 

PR'T(S)»0, 

CALI  SF  Ttr.(AMlNJ) 

!F(h«(xFN(>-X))JB,JT,2 

2 C'NTINUF 

A(l )aO.S 
»(2)«0.?92b«J2 
A(i)a)  ,7071  07 
A (0)30, 16866667 
6(1132, 

P(2lal  . 

R(5)b)  . 

«(U)s2, 

C(1)bO.S 

CtEFBfi.JRPBRJE 

0(5)31.707107 

CUlao.5 

DO  ] iFl.NDlM 
AU«fi,i)py(i) 

AUX(2,I)Pf'ERy(I) 

AlH(3«I)30t 
5 AIIX(6,1)30, 

IPECpO 

M«h*M 

I HLf  « • i 
ISTfp«0 

4 CONTINUE 

lFC(X4H»XtN0)*H)T»6,S 

5 CONTI NuE 

6 CONTXNUF 


Inlet 

I 7b 

inlet 

176 

INLET 

177 

INLET 

178 

INLET 

179 

inlet 

160 

INLET 

181 

INLET 

162 

inlet 

183 

INLET 

1 69 

inlet 

18b 

INLET 

186 

INLET 

167 

INLET 

166 

INLET 

169 

INLET 

190 

INLET 

191 

INLET 

192 

INLET 

19J 

INLET 

194 

INLET 

195 

I NLt  T 

196 

Inlet 

197 

INLE  T 

198 

INLET 

199 

INLET 

200 

INLET 

201 

INLET 

202 

INLET 

201 

INLET 

204 

Inlet 

205 

INLET 

206 

INLET 

207 

INLET 

206 

INLET 

209 

INLET 

210 

INLET 

211 

Inlet 

212 

INLET 

211 

Inlet 

214 

inlet 

215 

inlet 

216 

INLET 

217 

INLET 

216 

INLET 

219 

INLET 

220 

Inlet 

221 

inlet 

222 

INLET 

223 

Inlet 

229 

INLET 

225 

INLET 

226 

INLET 

227 

35 


BEST  AVAIU8LE  COPY 


»»*f NO.» 

I t~0« 1 
7 CONTINUE 

call  se»(H$.«) 

C‘LL  tPilHTE  f N I M|_t  T I K,nS.OInFLO.Y( AhInI.BNK.lT) 
IFLAGl*X/DfLTB 

IFUFLAGi.ht.tFLAGJ. and. 1TABLE. EO.l)  Call  table 
I»LaG2»TELaG| 

IF (PNmT(5) )A0.8.A0 
» CONTILUF 
I T t B T ■ 0 
9 CONTINUE 

ISTfPA JSIEPAl 
J»1 

to  CONTINUE 

A J* A ( J ) 

CJ*r( J) 

on  ii  i«i.nojh 

»l«M»OfPYf I) 

»?*»J*(MI-ajAAU»(».T)) 

Y(I)«V(I).B? 

It  AU>  ((,.  I )1A|J»(A.  I)*P?.CJA«1 
JMJ-A)12.i5il5 

12  continue 
J*Jt  I 

if l J-3) I J. |U. 13 

13  continue 

1«  continue 

call  SETEQ(ANINI) 

GO  to  in 

15  CONTINUE 

IE(TTESTJ1»«16.20 
I 6 CONTINUE 

00  IT  iAltNOI" 

17  AU*(U.I)1Y(I) 

IT1STA1 

ISIEP*lBTEp«ISTFP-2 
IB  CONTINUE 
I"LF*lNLE»| 

*■*»" 

H»0 ,S*H 

do  ip  isi»*ni» 

V ( I ) ■ V ( 1 • I ) 

0t«v(I)«4uif(2*X) 

I*  AU»(6.I)»Au«(3.T) 

GO  TO  9 

20  CONTINUE 
I»U0«ISTEP/2 

1F(  IBTtP«I"00«I*<00)2|  ,23*21 

21  CONTINUE 


INLET 

228 

iNLtT 

229 

INLET 

2 JO 

INLET 

2 J 1 

INLtT 

2J2 

Inlet 

233 

INLET 

2J« 

INLET 

235 

INLET 

236 

INLET 

237 

INLET 

238 

INLET 

239 

INLET 

240 

INLET 

24! 

INLET 

24? 

INLET 

243 

INLET 

244 

INLET 

245 

INLET 

246 

INLtT 

247 

INLET 

248 

inlet 

249 

INLET 

250 

inlet 

251 

INLET 

252 

inlet 

253 

Inlet 

254 

INLET 

255 

INLET 

256 

INLET 

257 

INLET 

258 

INLtT 

259 

INLET 

260 

INLE  T 

261 

INLET 

262 

INLET 

261 

INLtT 

264 

INLET 

265 

INLtT 

266 

INLtT 

267 

INLET 

268 

INLET 

269 

INLtT 

270 

INLtT 

27! 

INLET 

272 

Inlet 

273 

INLET 

274 

INLET 

275 

INLET 

276 

INLtT 

277 

INLET 

278 

INLtT 

279 

INLET 

280 

I 


36 


C‘Ll  Sf  TED(A"IM) 

Inlet 

26! 

DO  ?1  Ul.NOt" 

Inlet 

26? 

*u«(iin»y(n 

Inlet 

263 

22  *i'«  r 7.  i ) 

inlet 

209 

on  id  « 

INLET 

265 

21  continue 

INLET 

266 

OfLTsO, 

INLET 

267 

do  i*i,\ni" 

INLET 

266 

2d  DE  L t*DE 1 T » aij  X (8,  I )•  tgg(  lulllli  Ij»y(  I ) ) 

INLET 

289 

IMotLT-PWHT(U)  ) 28. 26,25 

INLET 

290 

25  CONTINUE 

INLET 

291 

IF(INLF-|0)2»..H,.l* 

INLE  t 

29? 

26  COntJNuF 

INLET 

293 

DO  ?y  1*1.601" 

INLET 

29a 

27  6iura>Z)«6ux(St>) 

INLET 

295 

ISlFP»lS’f P»ISTFP"6 

I "L  ' T 

296 

l«X.H 

INLE  T 

297 

IF  NO«0 

Inlet 

296 

GO  TU  IP 

INLET 

299 

26  CONTINUE 

Inlet 

300 

CALL  SFUO(*Ml«ir  J 

inlet 

3oi 

DO  ft  I*|. dOI- 

INLET 

302 

*'ix(i«n«*in 

INLET 

301 

*D«(2«l)»OfP*(I) 

INLET 

309 

*"*f  |.  MXAllKft,  T) 

INLET 

505 

*r  ii*»ux(s, i) 

Inlet 

306 

26  CE"t  ( I ) *ADX  ( 7.  I ) 

INLET 

307 

C»L1  Stilus. 1-H) 

Inlet 

308 

CALI  TP.HTt  (NINL£T,*..>,."S.UtNFLO,».  a"Int.RN6,NT1 

Inlet 

309 

I*i»r.t»rx-M)/nEi.T6 

inlet 

MO 

I'(TFLAGl.NF.lALAG2,AN0,ITA8LE.EG,n  CALL  TABLE 

INLET 

31! 

IFL*G2*IFlaS1 

INLET 

312 

IF (p»»7 [f  > JUO.  Vl.«0 

INLET 

311 

10  CONTINUE 

INLET 

319 

on  ii  i«i,M)i« 

inlet 

315 

* ( I l»»U« ( 1 . I ) 

INLET 

516 

11  PFHV ( I ) IAUI (2. I ) 

INLET 

317 

l«tr«l"LF 

INLET 

316 

IF(IE'-0)}2. 12.16 

INLET 

319 

12  continue 

inlet 

320 

I"LF«I"I  F-I 

Inlet 

321 

ISTFP6IS1FP/2 

inlet 

122 

inlet 

323 

I» ( tPLF16.ll,]} 

inlet 

329 

11  ContT^UF 

INLET 

325 

I"00*ISTEP/2 

Inlet 

326 

JE( TSTtP-I "OD»I "00)6.16 « 6 

INLET 

327 

16  CONTINUE 

INLET 

328 

IF(OtLT-0.02«P»"T(6) ) IS. 15 ,6 

INLET 

329 

15  CONTI"OF 

inlet 

330 

I"LF»I"lF*1 

Inlet 

331 

ISTfp»ISUP/2 

inlet 

532 

"•"♦h 

INLET 

333 

00  TO  6 

INLET 

359 

lb  CONTI NUE 

inlet 

355 

I"LF*l I 

INLET 

336 

CALL  StTEOtAMUn 

INLET 

337 

(.0  TO  16 

INLET 

336 

17  CONTINUE 

INLET 

339 

I "LF* 1 2 

inlet 

390 

GO  TO  ]« 

inlet 

391 

16  CONTINUE 

Inlet 

39? 

I"LF*I 1 

inlet 

393 

16  CONTINUE 

inlet 

399 

CALL  If  A ( "5  * X ) 

Inlet 

395 

CALL  TP.6Tf{Nl,(LfT,»,"S,alNrL0.»lA«lNl.6NK,NT) 

inlet 

396 

If  LA(,l««/Of  LT« 

Inlet 

397 

IF(TALAGl.Nf,IfLAG2,AND.ITA«LE.tO,l)  CALL  TABLE 

Inlet 

396 

IFLAG2*tf LAGI 

INLET 

399 

««  continue. 

INLET 

350 

PETupn 

inlet 

in 

(NO 

Inlet 

152 

37 


A 


best  available  copy 


S'.mnOUTTNf  SFTEH(**IN> 

INLET 

3S3 

C ROUllNf  TO  SfcUp  THf  EUuATlO^b  PflW  THfc  WIGHT  * A nn 

Sint  OF  Tut  tniuTIQMS 

I*ttT 

3b9 

C HOTJOV  (Ml  TO  Tut  BEN*  Of  Thf  Tt«“S 

1* 

ThE  tOUETION  OF  *07 J 0 

JNLtT 

3b5 

»fn  L»lE''G’**LTNtL«>NiN><.LF 

l NLt  7 

356 

COMl'Oi./l.'U«S/N!)f;.nlNLIT>ICM<J).tSl<i).(dH)L<7. 

7)>B(7.7).DC7,7)> 

I NUt  T 

3b7 

1 »(7.7).’7(7*7)»»l7.7).V(7>7)»n(7.7)."S,Mtt,Hr7,7),IC.IS.*»lM(J)> 

INLET 

3bft 

XH"lMK  1)  .LlN.QXf  3)  .QlMFLOt  ARAVfLf>GTM(3) 

I'-LtT 

3b9 

^nruu\/N7UMi/y(5)»nFOv(S)»X*‘<T»I»*TfZFTA*MH 

INLET 

360 

Ca'1HOK/DiUM?/0Kn«797)tDX(l«7tT)»HX(S«7v7)i 

i««(J.7,7).L*(J,7,7).N*(i 

INLET 

361 

1*7*7) 

INLET 

362 

CnM*(iN  /NIJM3/A0*  T*  ARVtBbTA 

INLET 

3&3 

f J*  a ) 

I *L  t 7 

369 

UI^F'jSIPn  amIM3) 

INLET 

365 

0«32.2 

INLET 

366 

DO  ??0  *il»j  . J 

INLET 

367 

nn  11*  r«t,rt 

INLET 

366 

\\9 

H9nrsi*j)*p, 

!Mt7 

369 

220 

COxTl^UF 

INLET 

370 

CALI  »fA(HS«X) 

I*Lf  T 

371 

hhihS 

INLET 

372 

c FINO  Thf  HAt  A 0 1 A 

INLET 

373 

HH«v  f ) 

I NLt  7 

37*4 

A«Av»A«V*(l#tRFTA*HR) 

INLET 

375 

0 r • n , 

I L fc  7 

37b 

C SfT 

l"»  Fii'JATlOKS  Fp»  fc*Cw  I*LtT 

INLt  7 

377 

pn  ino  ’Jl»i*MNi.tT 

I NLt  7 

3 7 A 

A“I»<(  Ml  *99999999999, 

INLE  t 

379 

I»Lt7 

360 

OTxrjT^OO 

INLE  7 

361 

ir«iCH(\n 

INLt  7 

362 

X $*  TSf (Ml ) 

INLET 

363 

L(  •*. 

INLE  7 

36« 

no  9S  l« i • is 

INLET 

365 

DO  9U  Jal'lC 

I NL  t T 

366 

M J#J)»*.XfNl#I#J) 

INLtT 

367 

l x(sl*l*j) 

I N L t 7 

366 

L F »i  f ♦tf l*J)/fFLOAT(XC)) 

I *•  L t 7 

369 

9u 

H(  1 • J ) a*x ( M • I « J ) 

INLET 

390 

9b 

CONTI^UF 

INLET 

39) 

CALL  LtVFL 

INLtT 

392 

A3*6. 

IME  7 

393 

i««A, 

INLET 

39a 

AF»n  . 

hut 

395 

on  9 7 Is l « IS 

INLtT 

396 

<1**0, 

INLET 

397 

Dl«0. 

INLET 

396 

00  9h  Js  1 • I C 

Inlet 

399 

I • J)/ (FLOAT rm*U) 

INLET 

900 

o(  I.  J)»nx(Ki* I* j)»Mf i,j) 

INLET 

901 

IF (DC I*J) .L’ ,0.)  0(T*J)»0«9M 

INLET 

902 

Inlet 

903 

IFUtl.JJ.lT.O.)  »(  T.  J ) ■ 0 , 0 0 1 

inlet 

909 

17(1, fO. 1)  iS«*8*»(TtJ) 

INLET 

90S 

38 


IF(I.EO.IS)  AH»ABTA(I,J) 

Inlet 

406 

«6 

A AAA  A* A f I • J! 

inlet 

907 

IF  ( AA.LT.AMIN(M)  ) *H!N(NIJ.AA 

INLtT 

906 

97 

»E»Af  »iJL/*A 

INLET 

909 

A“I»,i(ntj«amin(nI) 

INLtT 

910 

At«l ,/Af 

INLET 

911 

IFIThT.eu.I)  call  n t 1 

INLET 

912 

IF(i»T,Ft.?)  Call  «tj 

INLET 

911 

ifij.t.fu.sj  call  «tj 

INLET 

914 

DO  luO  1*1. IS 

INLET 

9J5 

o.i  i jr  j*i  . rc 

INLET 

916 

H«(Nt.I,J)»M(I,j) 

INLET 

917 

139 

«»(MI.I.J)»*(I.J) 

Inlet 

910 

190 

CONTINUE 

INLET 

9 1 9 

BNKfllI,?)*AE/(2,*LEl«(l,/(AR**2J.l,/(AS«*2))*OQ»0O 

inlet 

9 £0 

RNK(nI.J)*CAaE/LE*(hb.hS) 

inlet 

921 

Dll  B5  I*  1 t IS 

INLE  T 

922 

AC*n. 

INLET 

921 

DO  AH  J* 1 . IC 

INLET 

929 

09 

AC*AC»A(I.J) 

INLE  T 

925 

I'D  *J  Jil.IC 

INLET 

926 

03 

RNK(nI.iO*RNk(NT.R)*AF/<LE»aC>»G*N(I.J)**2*aBS<«CIfJ>*0G)» 

INLET 

927 

1»(I.Jl«i;G/(J.20A*0n,J)*»0.JJJJJ*AtI.J)A*?)*L(IlJ)»B(lFJ) 

INLET 

920 

05 

CONTINUF 

INLET 

929 

B N K ( N I , 1 )«.HN«(NI.J1.RNA(NI,J).HNK(N1>U1 

inlet 

910 

D(Rv(NI  )*«..K  (NJ  , 1 ) 

INLET 

911 

: MnO  Thf  RELATIVE  RANK  OF  TERMS,  NORmAlIIE  BY  T ME  LARGEST  TEH*. 

INLET 

932 

X"A*»0, 

INLtT 

931 

DO  101  I«1  ,9 

INLET 

939 

101 

1^  (AtlSCR^^C^I  • I ) ) *CT#X^AX)  XrtAX**BS(«^K(NI#  I)  ) 

INLET 

935 

DO  102  T • 1 * 9 

INLET 

936 

10  2 

T)«lO0#*RNK(*J*  J)/X«AX 

INLET 

937 

lo  n 

CONTl^ue 

INLET 

936 

DE»y(NlNLET*l)«QT/AR/kVt01NFL0/ABAV 

INLET 

939 

RETURN 

INLET 

940 

£*0 

INLET 

441 

SUHROI'T  rwf  TP-.RTE  rNTNLET*X*HS*OlNrLO»  Vf  A»lNT«RNK|  NT) 

INLET 

442 

subrouttnF  to  »r:te  hydraulic  Information  on  tapes 

INLET 

441 

OlMFNStnN  RNK(S.II)  ,V(S)«AHINI(  J) 

INLET 

444 

MOURS*«/J*lOO, 

INLtT 

445 

NT»nT*1 

INLE  T 

446 

00  100  Nl*l«NlNLET 

INLtT 

447 

IUNTT*NI*A 

INLET 

440 

»aT(»I)/afi«I(ni) 

INLET 

449 

too 

► RITE<  IliNin  hours,  mj.oInFloi  TIN  Inlet*  I ).V.r(Nl) , (Rnk(nI.J)  ,j«j ,u) 

INLET 

950 

return 

inlet 

951 

(NO 

inlet 

952 

39 


BEST  AVAILABLE  COPY 


SUBROUTINE  LtVfL 

iNLtT 

RSI 

c T*is  MGiiTInT  COm^'*TFS  LfVFLS  THROUGHOUT  Tm£  INLET  ASSUMING  LEVEL 

INLET 

uS« 

C ARfc 

l t nf xr  run*  rxy  to  sfx 

inlet 

RS5 

•*►*1  L.lE'-GTm.lIN.LX.N.NX 

INLET 

uS6 

COFMuN/^uMg/si.g.LUjLET.ICMli)  .IStn)T0RtL<7.7>.B(7i7).D(7,7). 

INLET 

«57 

1 x(7,7),N(7.7>.*<7,7).vC7.7).0(7.7).NS.Hl).H(7.7>.IC.tS.XMlM(J). 

Inlet 

«SH 

l B" I n I ti)  .LIN.Qxt  J)  .iJlNFLO.XgXV,  LENGTH!  J) 

Inlet 

aS9 

on  ?n  jpi.jc 

INLET 

4460 

XL«n, 

INLET 

R61 

DO  10  1=1.  IS 

inlet 

Ub2 

10 

XL*XL*L(I.J) 

inlet 

R65 

X»«L(1. .!)/?. 

inlet 

utu 

x|  1 ,J)»h5*(HB-hs)/XLRXX 

INLET 

4465 

DO  11  I»4.IS 

Inlet 

<466 

XX«(L(I.|.J)*L(I.J))/?,XXX 

inlet 

Ub7 

n 

Hll t J)*MS« (MB«M5) /XL*** 

inlet 

Ub» 

20 

CONTINUE 

INLET 

4469 

Rf  Turn 

INLET 

«70 

tNO 

INLET 

«71 

SUhWtlUT  TNF.  SFX(mS.TTmf) 

INLET 

«72 

c This  Subroutine  OETEN'iMfs  ThE  FORCING  SEX  level  EIThEs  FRO" 

inlet 

<473 

C EOUXL-I I"t.SEHlES  DXTX  { I F XVXlLXBLt)  OR  BT  SlNUSOOlXL  FORCING, 

inlet 

1*7  u 

COMMON  /niihS/XO.T  * AR « Bt  TX 

Inlet 

«75 

DIMENSION  r(52) 

Inlet 

«76 

UNPUN*  1 

inlet 

«77 

I‘  (Il’I.NF,  1 ) GO  TO  1 n 

INLET 

44  7* 

RExntN.I)  TDEL.NF>Ts 

INLET 

4»?9 

1 » UNMXT (XNX.F 0,2, 61. T J) 

inlet 

*460 

TDEL«Tl)FL*6(l. 

INLl  T 

«61 

C MF.XO  six  LEVEL  EGUXL  TINE  SEhIE*  OXTX  T«E  FIRST  TI*>E  SEX  IS  C»LLEO 

inlet 

4*b? 

c IF  NFMS  IS  r,«F*TER  thin  I 

inlet 

*4  rt  1 

lFl'ipTs.GT.l)  RFX0(S,2)  ( T ( J ) , J«  1 , NpTS ) 

inlet 

44  0 44 

i FONmjT ( «F 10. SI 

INIET 

44  b 5 

IFinpTS.GT.I)  «PITEI6.3)  (TtJ).JPl.NPTS) 

INLf  T 

4*66 

J FOR-xTdx,  !6F»,2) 

INLET 

44b  7 

MPNPTS*  1 

INLET 

44  6 * 

n?=npts*2 

INLtT 

4*69 

T(Nt)P»m 

INLET 

4490 

TfN«>«Vf2> 

INLET 

449  1 

io  ifimpts.lt, i>  go  to  ioo 

INLET 

4492 

C INTERROi .TE  in  time 

INLtT 

4493 

I T*  T I ME / T 

INLET 

44  V 44 

xt«ti«e-IT=t 

Inlet 

449S 

JxxT/TOEL 

INLtT 

4496 

J"J*  I 

INLET 

449  7 

HS*y( J)*( (T( J*1 )«T( J) )•(»!•( J-l )PTDFL)/TOfL  ) 

INLET 

«9* 

return 

INLtT 

4499 

C DETERMINE  level  if  SE‘  level  FLUCTUXTIUN  IS  SlNUSOOlXL 

INLET 

soo 

100  HS«»0=  SIN(2,XJ,l«lSg.TlME/T) 

INLET 

501 

RETURN 

INLET 

502 

END 

INLET 

503 
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SUBROUTINE  MELMTTheum.aH.CORL)  INLET 

C ESTIMATE  The  KUT.hu  mELmHULT?  PERIOD  INLET 

c of  tke  tut t/iu*  system  (neglect  friction)  inlet 

real  L.LEnGTm.lIN.LY.N.NA  I Nl  E T 

CnMMON/NUMS/NltG«NIllLETfICM(J).ISF(  I).UR.LC7.7)«Bt7>7)«D(7,7).  INLtT 

1 *(7.7),1.(7.7),«(7,7).V(7.7).U(7,7).MS,HB,M(7.7)#ICtIS.AMlNI(l),  INLET 

IP“InI  (1)  .LIN.OXfl)  .rJINELO.ARAY.LE  NGTMl J)  JNLE  T 

PlMfNSION  CORL(T)  INLET 

C USE  EIYE  ITERATIONS  TO  OBTAIN  TmE  ESTIMATE  INLtT 

DO  1000  lIll.S  INLET 

SUMiP,  INLET 

00  100  NMI.NlNI.tT  INLET 

A“1niA»ini (NN)  Inlet 

100  sum«sum*anIN/(lenGtm(nn)»CORL(nn)J  InLET 

TMtLMSi.RJ.UISR*  SORT(Am/G)/  SORT(SUM)  INLET 

C tSTIMATE  The  HELMHOLTZ  PERIOD  INLET 

no  toi  nnpi.nimet  Inlet 

c estimate  tme  inlet  length  correction  due  to  radiation  inlet 

101  C'>ML(NN)P-HMINI(NN)EJ,lUlS«PALOG(l,UlS»PeMINI(NN)/(  SORT(  INLET 

l3i.?*»MINI(NN)/RMlNI(NM )*ThELM) ) INLET 

1000  CONTINUE  INLtT 

C convert  tme  HELMHOLTZ  PERIOD  TO  HOURS  INLET 

THtLMPTMELM/JfcOO,  INltT 

return  Inlet 

END  Inlet 
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51 1 
517 
511 
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515 

516 
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5 27 
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525 

526 
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S'lMROUTINE  «TI  Inlet 
C This  SUMRUuTINE  "EIGHTS  Th(  FLU"  IN  t ACh  Section  SO  THAT  FRICTION  INLET 
C IN  THAT  SECTION  is  “INJMI7ED,  THIS  “fANS  THAT  AT  EACH  SECTION  FLO"  IS  INLET 
C ALLOnED  TO  REDISTHiruTe  ISTStLF  Throughout  THE  channels  TO  minimize  FR  INLtT 
C m0"E v E m » F Lf"  pEmPenoICiiLaR  TO  THE  CHANNELS  IS  ASSUMED  TO  BE  Small  *no  Inlet 
c flo«  is  not  included  in  the  equations  of  notion,  by  minimizing  fricti  inlet 
c routine  gives  an  upper  l imi t fur  ray  LtvtL  fluctuations  and  inlet  veto  inlet 


»e*l  l.i £ngtm.ltn.l*.n,n«  Inlet 

CO“MON/MU"5/Nl,r,,NlNlET.lCH(J),Ist(!).aR.L(T.T).B(7.T).0lT,T).  INLET 

1 *(7.7),nC7,7I,.(7,t).V(7.7).D(7.7).«S.mB.M(7,7),IC.IS.AMInI(J).  INLET 
!M“INI( 11 .L IN,0«(5) .nINFLO.ARAY.LENGTMt Jl  iNLtT 

DIMENSION  C(iD)  INI t T 

op  i no  iii, is  Inlet 

su"rio.  inlet 

DO  so  Jll.lC  INLET 

Cf J)*»( T.J)**Z*(D( I. J)i*.!ll)/  INLET 

1 (N(i.Ji«"z«c*(Ni)*ti*u(t.Ji*LU.J))  inlet 

50  SUMcpSuhCtCIJ)  INLtT 

DO  NO  Jll.lC  INLET 

60  «fI.J)iClJ)/SUMC  INLET 

loo  continue  inlet 

return  Inlet 

ENO  INLET 
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507 

508 
504 

550 
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BESIJVAI1ABIE  COPY 


SiiUMfiiiT  In*.  "T? 

INLET 

552 

c ROUTINE  TO  DtTfRHlNF  T„f  r.KIu  *f1GhT1«g  FUNCTION  ASSURING  THAT 

INLtT 

SSI 

C FLU*  IN  i GIVE  N CHANNEL  IS  ThF  s*Mt  ALONG  The  INTI«F  CHANNEL 

INLET 

SSu 

C KO 

. TS  DlSTRlhUTED  IN  [HAF'NtLS  '0  GIVI  A hJNImuh  TOT»L  FRICTION 

INLET 

sss 

C FRICTION  I>i  This  RnuTlaif  .ILL  RE  SLir.~TLT  hIG»{H  THAN  IN  .Tl  AND  THE 

INLET 

556 

C IN 

This  S v S T t m Is  CONSISTENT  pith  THt  i rj  u A T TONS  OF  *OTlON, 

INLtT 

557 

Rt*l  L •LtNfiTH#LXN*L  W»N*N* 

Inlet 

558 

COMMijN/Kiu*-s/f,I#r,,NlNLFTt  ICM(i)  *ISf(3).UR.l(7.7).8(7.7)  .0(7,7)  * 

Inlet 

559 

1 A (7.  M.M7.7)  ,*<7,7).V(7*7>fQ(7,V),MS,Mtt,M(7.7>#  IC»I3»AHJNI<1), 

INLtT 

560 

1HMINI ( S) .LJN.Olf 3) ,OINFLO#A§AV#LtNCTH( J) 

INLET 

561 

C(20) 

Inlet 

562 

S"MC«0, 

INLET 

563 

DO  100  T«1,IC 

INLET 

566 

C(I)«0# 

inlet 

565 

no  no  j* i • is 

INLET 

566 

50 

c f I>*C(  T )♦  (N(  J,  !)*v?tQX(Nl)M2'(H(  j*  l)*l(  Jt  I)  )✓ 

INLET 

567 

1 (A(J*I)**£*(  Of J* I) )••• J1331) ) 

INLET 

565 

C<1)»1  ,/CU) 

INLET 

569 

100 

sLiHr«suMC»c(  i ) 

INLET 

570 

00  70  J«1«IS 

INLET 

571 

on  #,0  la  1 • rc 

INLET 

5 72 

60 

*(J.I)»C(!)/SU*C 

INLtT 

573 

70 

continuf 

INLET 

576 

RETURN 

INLtT 

575 

END 

INLtT 

576 

S' '09C*< T T -T3 

INLtT 

577 

C 

THIS  ROUTINE  A gsljH t s That  OlSCHANGt  Is  ElIUALLT  distributed  THROUGHOUT 

INLET 

5 78 

C 

The 

I“LF'  GRID  S Y S T F h , IN  GtNtRAL  Tuts  »IU  NOT  BE 

TRUE  BFCAUSt  IT  IS 

inlet 

579 

c 

DIFFICULT  Tfl  AfCuRATtLr  ORA . This  TYPE  of  grid  ry  eye  AND  FLOH  distrub 

INLtT 

560 

c 

changes  «ITm  Tint  IN  HOST  I NLt  T s , ThJs  ROUTINE  IS 

USEFUL  i*  GIVING  a* 

INLET 

581 

c 

VELOCITIES  and  Hay  i.EVel  fluctuations. 

INI  tT 

582 

c 

GRIDS  MTh  DEPTHS  LT  0,n|  FOOT  aHC  ASSuHtO  TO  HAVE 

no  FLOW 

INLET 

563 

hF«l  L >1  E ngTh.l TN<l**n.n« 

INLET 

566 

COHm(in/nU"s/M  ,r,tNiw|.t  T i IChO).ISE(J)«urfL(?>T) 

.8(7.7) .0(7,7) . 

INLET 

565 

1 »(T,T),NfT,n,u(7,7).Vt7.T).a(T,7),HS.HB,H(TiT).ICtIS.AHlNI(J)t 

INLET 

566 

IHHINI  ( St  (I.  INtOA(i)iOINFL0AARAY«LtNGTH(J) 

INLET 

587 

00  ? I»l . Is 

INLtT 

568 

File 

INLET 

589 

CO  1 J«1«IC 

inlet 

590 

1 

IF(0(I.J).LT.«,M>  y ■ X"  1 a 

INLET 

591 

IFty.Lt.O,)  «RITE(6tlOO)  NialS 

INLET 

592 

100 

FURHAT(///,5x. | ERROR  ..  INLtT  has  DRIED  UP  AJ 

INDICATED  In  nTS (•/ 

INLET 

593 

) NF  a IINLtTR (a lua  1 SE C T ION.  | a 1 u a // / ) 

INLET 

596 

IFty.LE.Oa)  STOP 

INLtT 

595 

00  S J«l a ic 

INLET 

596 

.( I a J ) ■ 1 a/y 

INLET 

597 

3 

IF(D(I,J),LT.O.nl)  H(I.J)PO, 

INLET 

598 

2 

CONTINUE 

INLET 

599 

be  Turn 

INLET 

600 

INI) 

INLET 

601 

SUhROUTT*^  TaRlF 

I*lfcT 

60? 

; KOUTInL  til  »MTt  4 T4f<LF  nt  IkST4M4lE  OuS  MVUR4ULICS 

INLtT 

603 

3t*L  L.l  t‘‘CT"*LT>  *l».n*nX 

iNLtt 

604 

CP""UN/hg*'s/Nl  1»  LFT*  1C**(  J)  * ISE  ( S)  *«B*L(7t  7)  *»<7*  7)  *0(7*7)  * 

INLtT 

605 

1 A(7«7).^f7*7)»w(7.7)*V(7#7l»0(7#7).HS,wb»H(7.7)*IC*IS»*RlNl(3)* 

INLtT 

606 

JMLlN.ilXf  J),OlNELO**<UY*LtNGTM(  J) 

I nl  t T 

607 

CONMtlN/M.1**  1 / Y (5  1 *OlfiY(S)  • X • NT  « 1 NT  * / F T 4 * MH 

I NCI  T 

606 

CO"HnN/Nu»2/Rx(X*7,7)*OX(3*T*7)*HX(J*7>7)*Hl(}i7*7)*LX(l*T|7)*NX(S 

IMIT 

609 

t i 7*7) 

INLl  T 

610 

Cm'Mp.i/Mlittu/RMl  (J#  4) 

INLET 

611 

DI^FnSIDn 

iNLf  T 

61? 

li4T4  N4ME/»hv(FpS)  ,*,H4(FT?)  .*H„EIGMT  ..HLEVtL 

Inlet 

613 

HHiiax/ihOO  # 

INLET 

614 

»«i7t(».n  h“s 

INLET 

615 

1 

INLtt 

616 

is«*  tri«F.  hou«s  • i.fe.j) 

Inlet 

617 

CO  100  MU  iMM.tl 

INLtt 

61« 

* «Ite  <l>.  IP)  M .xs.HU.y  (M  ) 

I NLE  T 

619 

10 

tOB-»T  (/,iox.  r I MLt  T (•  Tlt/t  lnx.  (St*  LtvEL*FT»[.F7,2./.lOX,(g4Y 

LEVE 

iM.fcT 

6?0 

ll.t7»l,t7.i./.10X.(nISCM»“GF.CFS»(,fl0.o«/,ai,  [CM4NNtL 

SECT 

INLtt 

6?  1 

lio  1 ? S u S 6 1 ) 

I NLt  T 

6?? 

ir«  tcM(M  > 

iMt  T 

6?3 

JS«TSfc (Ml ) 

I NL  t T 

6?u 

Do  n J«  i • I C 

INLtt 

6?  S 

D»*  3 I » i • T S 

INLtt 

6?6 

4<I.J)iNx(,.I,I,.i)»(0x(*.I.I.j)*NX(NltI*J)>»Mi(Nl.I*J)»*B8(Hx(Nl 

*I.J 

IMIT 

6?7 

1 ) )/f  71  1 »•»  l (•*  T ( [CM 

Inlet 

6?6 

I N L t T 

6?  9 

V ( I • J)«Y(>  Il<»xl  I.T.J)/4<I,J) 

INLE  T 

630 

3 

IT  ( »( I..I)  *Ll  ,n.M  ) v(  I,  J)*0, 

INLE.T 

631 

I>(J,tl.,n  •HITFC6.S0)  J.N»“t(U).(HX(BI,ItJ)>I«l.IS) 

Inlet 

63? 

•RlTfc (6.69) 

inlet 

633 

69 

»OB“»T ( /) 

IMF  T 

634 

*»]Tf  (s.SO)  J.K»hF(1).<¥  <!•.!}•  1*1.  IS) 

INLtt 

63S 

so 

(r(.«»T(ox,I?,lx,»s,?x,t.F  10,2) 

INLtt 

636 

-»i7M»,so)  (?) , (4(i i j) . in , is) 

INLET 

637 

• wlTF  (».S0)  JfS4Ht{Y)f[wx(h>T*ItJ)«l4)*IS) 

I NLt  T 

6 JB 

4 

CC'NTINup 

INLtt 

639 

•»lTt(b*S.)  (HN«(NJ,ir)*II*U«) 

INLtt 

640 

*9 

F0B«4T(SX,  (TfcNP  4CC«I«»7,l't  CONV  4CC»I.f7,tl(  H£  40»  ( • 1 7 , l . I EMC* 

I NLt  T 

641 

1 1*77,1) 

iNLtT 

64? 

Yf>AO«V  (*'!)/ A-lNT(M) 

INLET 

643 

.»1tf(b,»M  yH40  *4-I*,I(Nl) 

INLET 

644 

6 1 

/Platts!.  |“14N  VELOCITY  47  THE  HIHJHUH  4»t4  SECTION*  (*77*2*  l 

FT/S 

INLET 

645 

llC  1.  1 4 •* I n*  (*F«.i,  1 ru  1) 

INLtt 

646 

100 

CONTINUE 

INLET 

647 

HE  t(i»N 

I NLt  T 

646 

t'O 

inlet 

649 
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SLiHBOUTTNf  CRIT(nT,OFLT« IUNJ 1,1 ,NCrCLtS) 

Inlet 

6b0 

c 

Subroutine  CHIT  COMPARE*  1 COHSFCUTIVF  function  point* 

INLET 

651 

c 

ANO  -RITES  "IODLE  POINT  IF  JT  IS  A CRITICAL  POINT 

INLET 

6b2 

c 

I "Lt  T 

653 

OlOfNSION  F(1«S)  >"ABX(S>  «TER"(a> 

lNltT 

bbu 

data  "araa/im  /,  "APHR/IH*/ 

Inlet 

69b 

Rt-TNl)  TU-IT 

Inle  T 

656 

NLlNESafl 

INLET 

6b7 

TFar/Spoo. 

INLET 

658 

«“m(».  im«i 

INLET 

6b9 

ro  i Nit,? 

INLET 

660 

1 

HF»0( IUNIT)  «. [FIN, J) ,jn ,bt • (T(RN( II) , Jill ,«) 

inlet 

661 

00  too  MS, NT 

INLET 

662 

HE  *0(1  UNIT)  X,(F(3,.t),J«llS),(TFHM(II),TIa|,a) 

Inlet 

663 

lF(a,LT.p|,nt*|A)  Co  TO  101 

INLET 

66  u 

lOUTlO 

W E T 

66S 

00  2020  !»•!.* 

INLET 

666 

• 4H« (To  • “AM  i 

INLET 

667 

IF  (F(2.IA)  • Ftl.IA))  2012,  2020.  201a 

INLET 

666 

201? 

IF  (F(3.I»)  - »f?.IA)l  2020,  20t<>.  20tb 

inlet 

669 

20H 

I*  (F  ( 3,  I A j . F(2.II))  201S,  20  IS.  2020 

inlet 

670 

C Critical  point  value  Found 

Inlet 

671 

2019 

I OUT  1 1 

inlet 

672 

"AHX(lA)  ■ "AR«H 

Inlet 

673 

IF(  11,10.1, AND, F(?,TA),GT,U.)  mSm.F (2, [A) 

INLET 

67<i 

I»  ttl.fO.I  .*NP,F(?,TA),CT.O.J  TIpX 

Inlet 

675 

I»(IA,£l3.1.AN0,F(?,TA).Lt,0.)  HSLlF(i.lA) 

INLET 

676 

IF ( T A, to, 1 ,4ND,F (2, IAT , Lt ,0, ) T2ax 

inlet 

677 

IF(TA,tn,3,AND,F(3,TA),CT,0,)  HBmF(J.Ia) 

Inlet 

676 

I»( IA.ln.5.ANO,F  (l,IA).CT.O.)  T Six 

Inlet 

679 

IF(ia.Eu.5,an0.f(3,ia3 ,Lt,0,)  "HL»F(4»Ia> 

Inlet 

6H0 

IF(TA.£0,3,ANr),F(3,TA).Lt.O.)  TU1X 

Inle  t 

661 

lF(Ti,in,«,ANn,r(?,TAj,LT,0,)  vfif{?.Ia) 

Inlet 

66? 

IF(  TA.f  Q.R.AND.r  (2,  IA)  ,CT,0.)  VFlF(?,  U) 

inlet 

663 

2020 

CONTINUE 

inle  t 

66<J 

UCr  ?02S  I A • t , s 

Inlet 

665 

F(.1,I»)  » F ( 2 * I A ) 

inle  t 

666 

2029 

F ( 2 , I A ) 1 F ( 3 , I A ) 

inlet 

667 

IF  (IOuT.ro, 0)  r.o  Tfl  too 

Inlet 

668 

C 

IFIx.LT,  (NC*CUs-2)»TF)  CO  TO  100 

INLE  T 

669 

nlInF.SinlI»ES»| 

inlet 

690 

IF(NLlNfS,OT. IS"1  CO  TO  100 

inlet 

691 

-RITE  (A  ,2101)  x.(F(i,!a),mark(ia),iaii,$) 

INLET 

69? 

100 

continue 

INLET 

693 

101 

HINT 

INLE  T 

69u 

A"Pmi«bh/“Sh 

INLET 

69b 

a"Plp"»I  /h*L 

INLE  T 

696 

F“"i  AHS(T3.t 1 )l3S0./TF 

INLET 

697 

P»L«  AHS(Ta.T2)Ai60./TF 

INLET 

690 

»»ITf (0, 1 01 1 ) A“«",PMM, VF , AMPL.PML, Vt 

INLE  T 

699 

-RITE (p. 1 1 t 1 ) TF 

INLET 

706 

till 

FUH“AT(  5x,  I tf ■ [ , F 7 , 2 ) 

Inlet 

701 

RETiibn 

I"LE  T 

702 

2101 

FUHxaT  (2FH,3,At,-3PFA.*>‘l,2(0PF7,3,Al), 

Inle  t 

703 

3PF  0,3.  At,  2 (FA, "3,  A|)) 

inlet 

70« 

looo 

F0N“AT(aX,a“TI"F,Sx,2«HS«aA,oHiNF  L0« • SX , 2H„r, 

Inlet 

705 

1 bx.3~vtL,Tx,iu(.,/,Sx,iHHKS,Sx,2MFT,Sx,aHACFS, 

Inlet 

706 

I OX,2HFT,i,X,3HFPS.NX.aMACF5,/) 

Inlet 

707 

1011 

F0N“AT(//,1X,  Ia  CRITICAL  POINT  VALUE  !,///« 1*1. 

Inlet 

706 

1 («AVl  PROPAGATION  (,/,  1 SX  . [All/AO  l.sx,  (Phase  LAC(OEC)  "AX  VEL  (* 

iNLtT 

709 

1//«2X,[hICm  "‘TER I.2X.3F 10. a,/, 

Inlet 

710 

1 2*,[L0H  »aTe«  t,2x ■ ]f 1 0,a) 

INLtT 

711 

END 

INLET 

712 
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SUBROUTINE  REAOIN  (X,V.YFAC,XFAC.XO.XF,INDC.»K.LN.IUNn) 

INLET 

715 

c 

SUbPGUT I NE  TO  RFAD  SOLUTION  TABULATION  FROM  FILE 

inlet 

710 

c 

iNLtT 

715 

DIVISION  v(*),  VFaCI*) 

I NLt  T 

716 

0T5i,5»l ,/fc n. 

INLET 

717 

RF  An  (IlINlT)  I|  V 

Inlet 

7 1 fi 

:f(k,lt,.i.i«io)  kk.2 

INLET 

7 1 * 

i«or  ■ o 

INLET 

720 

IF  (K*  » 1)  10,  10,  50 

INLET 

721 

10 

IF  rxo  . X . 0T5)  20,  50,  So 

inlet 

722 

20 

IF  IM  • XF  . 0T5>  SO.  25.  25 

Inlet 

723 

25 

K«  ■ 2 

INLET 

72« 

GO  TO  50 

INLET 

725 

10 

INDc  ■ 1 

INLET 

72* 

X ■ xF»c»tx  • XOl 

INLET 

727 

»(LN)  ■ rF4C(LN)»T(LN) 

Inlet 

726 

50 

RETURN 

Inlet 

72* 

F no 

inlet 

730 

SUBROUTINE  GHPmC<alaBL1>»L*PL2.0ELT,1UNIT.NI) 

INLET 

Til 

c 

INLET 

732 

c 

Subroutine  GRPmc  writes  PLOTTER  TXPE  FOR  GRAPHICAL 

INLET 

733 

c 

OUTPUT  OF  SOLUTION 

Inlet 

73* 

c 

INLET 

735 

DIMENSION  RL ( 2 ) , ISVM(5) 

Inlet 

7 36 

nt”FNSiPN  VLARLI.  ( 3)  ,»l EON (3, 6)  , ALABLI  <«)  « AL4BL2<«>  «SYP< 3)  «V(*J  |YM 

inlet 

737 

lCf*l.XXl20oO),v»I2onol.TT(*,2) 

inlet 

73* 

0*T»  TLXHLL/IOnnEIRmTS.  V , 1 ()MElOC  I T IES«  ,«M.r  T , 

*♦>3/ 

Inlet 

7 39 

OX  T 4 XLFGN/1 ONFLOW  (KCFS.10H)  ,SN  .tOMINLET 

vflo,io»city 

Inlet 

7*0 

1 (FT/S.3mic>.10»WAy  LEVEL(iIUMFT)  . 3H 

.lOHINFLOR  tlOH 

inlet 

7*1 

? .1"  ,1l)Noef4N  LEVF.lONLIFTl  ,1" 

.IOHLEGENO  t l OH 

INLET 

7*2 

3 . 5m  / 

Inlet 

7*3 

0 X T 4 HL/IOMOOSERVEO  R.lOHAY  TIDE  ! 

Inlet 

7«* 

0 4 T 4 ISVN/5.P.3.2.  I / 

inlet 

7*5 

0 XT  4 Tl(*,ll/10NTfMP0R»L  4/ 

INLET 

7*6 

04 1 4 TT(n,2)/10hCCIL  / 

INLtT 

7*7 

0414  TT ( 7. 1 ) /I OMCnNvrCTlVE/ 

iNLf  T 

7*8 

04l»  TT(7.2)/10H  4CC  / 

inlet 

7*9 

0 4 T 4 TT  (0.  l ) /I  ONPRE  SSIIPE  H/ 

INLE  T 

7*0 

DxT  A TT(6.2)/lO"HO  / 

INLET 

7*1 

l)4»4  TT(R.  1 ) /tONOOTTON  5TH/ 

Inlet 

752 

0414  TlfV,21/lo«ESS  / 

INLET 

753 

c 

Inlet 

75* 

c 

R£ x o infoRnaTion  to  oipfct  plotting 

INLtT 

755 

c 

INLET 

756 

c 

FIRST  f A R D 

INLtT 

757 

c 

XI)  - STARTING  TINf  OF  PLOT  (MRS) 

INLET 

758 

c 

XF  . ending  T]Mt  OF  PLDT  (NNS) 

INLtT 

759 

c 

SCaL x - Time  axis  SC ‘OE  IN  HOURS  PER  INCH 

INLET 

760 

c 

TLfl  - «INIH(JM  value  OF  TIO*L  heights  (FT) 

Inlet 

701 

c 

YL  - OVERALL  HEIGHT  OF  PLOT  (InLHES) 

inlet 

702 

c 

YL  SC  XL  • SCALE  OF  T104L  HEIGHTS  (FT/incm) 

INLtT 

7©3 

c 

YRO  • mINImUM  value  of  FlO"S  (THOuSanOS  of  cubic 

feet  per 

SECOND) 

inlet 

70* 

c 

YRSC4L  » SCALE  OF  FLOW  ( THOUSANDS  OF  CUBIC  FEET 

PER  SECONO/lNCH) 

INLET 

70S 

c 

INLET 

7 00 

c 

CAPO  2 

INLET 

707 

c 

VV(1  • minimum  VELOCITY  (FT/SEC) 

INLET 

7o8 

c 

VVSCaL  - SCALE  of  VELOCITY  (FEET  PER  SECOND/INCH) 

INLtT 

709 

c 

SCAL  ■ SCALE  FACTOR  FOP  total  PLOT  SI2E 

INLET 

770 

c 

10  • Not  EQUAL  TO  IERo  for  4 PLOT  of  inlet  DISCHARGE 

INLET 

T71 

c 

INLET 

TT2 
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IMxl.fl).)) 

IXLFT 

7/3 

ul*n  ( 5.200  1 ) «0.»2.SC*L»»*uO.*l.»L»C»i  .y«0.THSC»L.TVO.YVSC»Ll 

IXLFT 

77a 

l scaiF.m 

lXlfT 

775 

2001 

K)H*>»T(»>  l(l.5./.Jfl».S.IlO) 

IXLFT 

776 

»81Tt  ((..2002)  »0.l*.SC»L«.*tO.»t.»LSC»L.»«0.yfl5C»l..*VU.»V5C»lt 

IXLt  T 

777 

I SC»LF.I« 

IXLfc  T 

7 76 

2002 

F0**»T(///,5*.  (PLOT  IxFOh»»TIOX (•/ 

IXLFT 

7/9 

1 t«.«F  10,5,/.U.5Fin,5«I101 

IXLtT 

760 

c oeTfn»int  s*"h»l  sp»cixr. 

IXLFT 

761 

tIXTya«.25»SC»L«/(fi€LT/J»00,) 

IXLt  y 

76? 

•al r t ( o . 1 2 l 5 ) lIXTtp 

ixlft 

763 

1215 

F(1Hm»T(|»,  (LIXTYPa  (.  Ik) 

IXLET 

764 

c 

IXLFT 

765 

c 

PLOT  Itr.fXD 

IXLl  T 

766 

c 

IXLtT 

787 

C*Lt  StxBOl ( 1 . ..*L/2,..e. .2A.6«LtCtxu.o..») 

IXLfT 

768 

no  ?n  LX  ■ t>  5 

IXUT 

769 

jxD»  • n 

IXLtT 

790 

yPa.yL/2.».«-Lxy,2 

IXLFT 

791 

LLx«lSy“tl  x) 

IXLFT 

792 

C»ll  Syxho(.(o..yP*,o»,  ,l«.tix,o..»l ) 

IXLFT 

793 

S»"fl>  ■ »l.f6x(|,Lx) 

IXLFT 

79a 

S Y 1 2 1 • »lltX(p,Lx) 

IXLFT 

795 

sytji*  »Lt6x(j,ix) 

IXLFT 

796 

C»U.  Syx(jui.t  .u.ya.A.i  .sy“.B,.21) 

1NLFT 

797 

20 

coxrrut 

IXLFT 

798 

C PLOT  TlTLt 

IXLFT 

799 

C»Ll  Sy"6nt(i.5.-*u/2.-t...21.»L»HL1.0,,J2) 

IXLFT 

60  0 

C*Ll  Sy-XOlI  1.5,-yL/2.*l  .<*.,21  « *1  »B12  • 0 . . 12 ) 

IxLET 

601 

C Plot  »»fs 

IXLFT 

602 

YLOy.TL/2.y*LSC»l 

IXLFT 

Mo  3 

C»LL  »»TStO..-yi,/2..  IfcHVtLUCITy.  FT/S1C.16.YL.P0.YYVO 

IXLFT 

60  a 

1 .yy*C»l) 

IXLFT 

605 

C*IL  *<TS(..8.-YL/2..1  |M"EIG"TS.  Ft.  1 1 . VL  « PO • .VLO*  YLSC AL ) 

IxLFT 

606 

C*Ll  AxTSt0.eYL/2..a«Tl«F(  "«S.-P. («F.*0>/SCaI*.0..0.iSCAl«) 

IXLFT 

807 

I* ( tn.xf  ,0) 

IXLFT 

608 

'CALL  *«1SU«F-«0)/SC»I.»Y**L/2..10HH.O«,  *CFSyx)0»YL  .50.eYL/*.»YP 

IXLFT 

809 

ISC»L.*«SC»1  ) 

IXLET 

610 

ifitc,. io.o)  c*u  flatk  •F*«o)/scFL*.*yL/2..D 

IXLFT 

Mil 

IF(Tc.fo.O)  C»ll  FtOTf («F-»0)/sC»L*,yL/2,.2) 

IxLFT 

81? 

CALL  pL0Tl(«F-«0)/SCFL«Yyi./2.ti) 

Ixlft 

813 

CALL  »LOT(0..yL/2.,2) 

IXLFT 

6 1 a 

*►*0(1)  • l./YLSC*L 

IXLFT 

815 

y F»C(2)  • (i.OO|/y»st*L 

IXLET 

816 

yF*cU)  * yffcii) 

IXLFT 

617 

*F»r(«)  • i./yvsc»i 

IXLFT 

818 

y“tt5)  • y»»C(2) 

IXLfT 

819 

no  i2io  llox.o 

ixlft 

620 

12)« 

*nr(U)».ooj 

ixlft 

621 

»F»r  • i./sc»l» 

Ixlft 

822 

•0 

on  os  i ■ i»  * 

Ixlft 

823 

C I» 

lliau  DO  xut  PLOT  OIsCxarSF 

IXLFT 

62a 

IFMD.tO.O.oxO.I.eo.S)  00  TO  05 

IXLET 

625 

COHaYL/2.* ( 1.5)10 ,A 

IXLFT 

626 

C»U  PLOT  (0,,  0,.  5) 

ixlft 

627 

0*  a 1 

IXLfT 

626 

ISuaaO 

ixlft 

629 

Hl«IxO  TuxjT 

ixlft 

630 
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INUX  ■ 0 

IMLLT 

8il 

os 

CALI  RiAOlN  (X.Y.Yf  AC*X*  AC*xO*XF*I*»r>C*KK*X,IU*lT) 

IXll  T 

812 

GO  TO  (70*  80).  KK 

IMLLT 

631 

TO 

ir  l TnPC.C*  ,0)  Go  TO  t>S 

IMLtT 

8i« 

72 

I Sun* I Suri* 1 

1-Lf  T 

855 

1M  iSUh.GP.  J99fl)  I S('8«  1 998 

IMLI  T 

836 

XX ( TSUMTCX 

IMLtT 

817 

vv(  TSi|«)«v(!) 

1 FML  f T 

818 

lMl.GT.b)  TY(  I.SOH)»VV(ISUtt)fC0R 

!*Lt  T 

819 

P ( TSU8.iL. 1998)  Go  TO  80 

I MLt  Y 

6i»0 

GO  TO  65 

IMLLT 

64*1 

60 

XX ( ISUH4 l )«0 , 

IMLLT 

8«2 

XX{  ISU»42>»1 ,0 

I '-LIT 

84*1 

YV  USUH4l)*0, 

IMLtT 

8w4* 

W(TSUh42)«!  , 

iNLtT 

8u5 

c plot  cupvts  ( no  not  pl"t  if  ewuai  to  zt«o  throughout) 

IMLLT 

8*** 

J*(vyUsuh.<»)  .ro.o.a.ANO. 

IMLLT 

8**7 

1 YvriSon.ij.tO.n.o.AMO.rrUsuHI.tO.O.O)  to  TO  eS 

IMLLT 

8**8 

IF(I,GT.5>  00  TO  Km 

INLLT 

6**9 

C»ll  LIMt(xX.YV.ISuft.l.LIMTvP.I) 

IMLCT 

850 

to  to  ai 

ILL  t T 

851 

885 

C*ll  L I»'t  (xx.  TY,  isiio,  i < 0.0) 

IMLLT 

852 

C»Ll  PLOW  (XFmxO)/SC»LX.COK,J) 

IMLLT 

851 

C *LL  PLOH(I,,COB.21 

I MLt  T 

854* 

SY^M  )«TT  ( I . 1 ) 

iLLtT 

855 

5»"«*>»TTt  W> 

I ML  L Y 

856 

C‘lL  SY»»Ol{.2,2,CO».0.1*SYM,0,,20) 

IMLLT 

857 

8b 

CO.-.T  !'<UF 

IMLLT 

858 

C READ  PHD T 0 f vpt  HAV  Tinf  (OAtA  STARTS  at  bEGINRlRG  or  PLOT*SAHfc  HaTuR) 

IMUT 

859 

IHM.MA.t)  tO  TO  2010 

IMLLT 

860 

MF.An<S.l)  TOft.NpTs 

IMLLT 

861 

l 

»'0M»aT( J«*,F6.2,b», 1 J)  * 

iMLtT 

862 

1HmpTS.it. 2)  GO  Tn  2m« 

IMLI  T 

863 

IF  C‘jpT  S.tT  , 1 ) (YY(J),J«1,MPTS  ) 

IMLLT 

664* 

2 

M(H  M A T ( MF  Id,  S) 

IMUT 

865 

XX (mpts* 1 )«0  , 

IMLLT 

866 

XX(MPTS*2)*1 . 

Imlit 

867 

T»txPTS*l)»0. 

IMLLT 

868 

yv(ppT$*2)*t . 

IMLLT 

869 

DO  1 J»l .MPTS 

IMLLT 

670 

YT(  J)*YY(J)»YFxCai 

IMLLT 

871 

3 

XX(j)»(TOtL/hO.)»«FAC»(J*l) 

IMLLT 

872 

C*LL  PLOtfxX(l>,YY(l),J) 

IMLCT 

873 

C*LL  LIMKXX.Yy.MPTS.I.O.n) 

IMLLT 

8 74* 

C»LL  plot(xx(Mpts/2) .yyCmptv?) . J) 

IMLLT 

875 

call  plot (xx(npts/2) .yy(NPTs/2)*,Ts,2) 

IMLLT 

676 

C»LL  SYMBOL ( XX (NPTS/2J*. ).YY(NpTS/2> ♦.?*,, I,»l,0..1?) 

IMLLT 

677 

20)9 

CALL  PL'’T((XfMXO)/SC»L**,'.»0..«S> 

IMLCT 

676 

«C  t(IPN 

IMLLT 

879 

LMO 

IMLLT 

660 
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